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Studies of past human-landscape interactions rely upon the integration of archaeological, biological and 
geological information within their geographical context. However, detecting the often ephemeral traces of 
human activities at a landscape scale remains difficult with conventional archaeological field survey. 
Geophysical methods offer a solution by bridging the gap between point finds and the surrounding 
landscape, but these surveys often solely target archaeological features. Here we show how simultaneous 
mapping of multiple physical soil properties with a high resolution multi-receiver electromagnetic 
induction (EMI) survey permits a reconstruction of the three-dimensional layout and pedological setting of 
a medieval reclaimed landscape in Flanders (Belgium). Combined with limited and directed excavations, the 
results offer a unique insight into the way such marginal landscapes were reclaimed and occupied during the 
Middle Ages. This approach provides a robust foundation for unravelling complex historical landscapes and 
will enhance our understanding of past human-landscape interactions. 



I n landscape archaeology, geophysical methods are increasingly being applied to conduct large area surveys 1 . 
I Although these techniques allow obtaining high-resolution archaeological information at a landscape-scale 2-4 , 
I they often neglect natural landscape variations. To fully understand the driving mechanisms behind human 
land-use, the integration of pedological and geomorphological information in these prospection stages is crucial. 

In Europe, one of the most characteristic examples of past human-landscape interaction, is the reclamation of 
wetlands and forest that followed the urbanisation of the historical County of Flanders (Fig. lb). During Medieval 
periods, these lowlands experienced a considerable population growth, making the County one of the most 
densely populated areas in Europe between the 11 th and the 15 th centuries 5 . To meet the demands of emerging 
cities, such as Ghent, Bruges and Ypres 6 , entire natural landscapes were reclaimed, transforming these into a 
landscape of dynamic exploitation 7 . The ruling Counts of Flanders spearheaded this evolution by endowing 
feudal lords with lands, while abbeys were deliberately installed on marginal land all over Flanders. 

More recent, natural processes (e.g. flooding) and modernisation have altered the cultural landscape of 
medieval Flanders, leaving the nature and exact layout of these designed landscapes largely unknown 8 . 
Although social formation and economic exploitation of reclaimed landscapes during the High Middle Ages 
is, to some extent, documented in historical and cartographical sources, archaeological evidence remains scarce. 
One of the main reasons for this scarcity is that traditional investigations have focused on recovering architectural 
remains rather than situating the structures within their broader environmental contexts. By focusing on elite 
residences and monastic buildings, archaeological investigations have contributed to the biased image of medi- 
eval settlement landscapes that still prevails. 

The situation was no different for the Cistercian abbey of Boudelo located in the north of the County of 
Flanders (Fig. lb). In 1 197, a small religious community settled in this area 9 , which at that time was an outback of 
the County, dominated by marshes and wetlands. Early historical accounts, such as those by monks from the 
Abbey of Clairvaux in Bourgogne (France), give testament to the harsh environmental conditions. They describe 
life at Boudelo as pauper and misserimus: poor and full of misery 10 . However, the community improved its 
property by cultivating the surrounding land. Between the 13 th and 16 th centuries, the monastic estate expanded 
to over 1000 ha, making the community one of the leading cultivators in the County. After religious as well as 
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Figure 1 | Location and historical context of the site, (a), The location 
and elevation 32 of the study area with indication of the excavated Abbey of 
Boudelo (1). The dashed line represents the boundary between the 
palaeolake (to its south) and the coversand ridge (to its north), 
(b), Location of the studied area (red dot) within the late medieval County 
of Flanders (orange boundary) and the current state boundary of Belgium. 
Maps composed in ESRI ArcMAP 9.3.1. 

military struggles and successive floods, the monks were forced to 
abandon their grounds in Boudelo in 1578. They found refuge in 
Ghent and the monastic buildings were sold and dismantled. 

Although the extensive reclamation strategy of the abbey is 
attested in historical records mentioning embankment, drainage, 
stockbreeding, extraction of peat and clay, and the production of 
building ceramics, it remains hidden in the archaeological record. 
Excavations in the 1970s revealed the remnants of the cloister range 
and the abbey church, leading archaeologists to interpret the abbey 
grounds as limited to the coversand ridge 11 (Fig. la). The border of 
the monastic precinct was believed to coincide with the edge of a Late 
Glacial palaeolake, which was, and remains, a waterlogged envir- 
onment dominated by peat and lime-rich lacustrine deposits 
(Fig. 1, Supplementary Fig. 1-2). However, in 2011 an electromag- 
netic induction (EMI) survey of this wetland area drastically altered 
this interpretation, as it unveiled traces of the abbey's outer court that 
once was part of the monastic precinct; a previously unknown 
designed landscape from which the monks directed their cultivation 
of the surrounding area. 

The multi- receiver EMI sensor we used, allows simultaneously 
recording the apparent electrical conductivity (a a ) and apparent 
magnetic susceptibility (kJ of four different soil volumes 12 . Where- 
as in non- saline environments a a can be directly related to soil tex- 
ture (clay, silt, sand), and is influenced by soil organic matter and 
water content 1314 , K a anomalies are often the result from ferrimag- 
netic soil enrichment 1516 , the anthropogenic disturbance 17 of top soils 
(e.g. ditches), or the heating of soil 18 (e.g. hearths, bricks). This con- 
trasts to other geophysical survey techniques, such as magnetome- 
try 19 and ground penetrating radar 20 , that each target only one 



specific variable and are inadequate to appreciate the full range of 
soil textural variability of surveyed sites. However, integrating elec- 
trical and magnetic soil properties in a single EMI survey allows 
targeting the pedological, geomorphological and archaeological var- 
iations jointly. By simultaneously mapping multiple soil volumes 
with a multi-receiver instrument, vertical discrimination potential 
is added to the electrical and magnetic data layers. 

Results 

The a a data revealed a complex environment with alternations of 
highly conductive lake deposits and peat accumulations (Fig. 2a). 
Apart from scarce sandy outcrops the measurements indicate a bur- 
ied marshland influenced by long lasting, saturated conditions. 
Several ditches were revealed as highly conductive anomalies sur- 
rounding the entire survey area and enclosing two smaller zones. The 
high conductivity of the ditches is caused by high concentrations of 
fine-grained sediments (clay) and peat that contrast with the more 
resistive underlying sand. The extent and layout of the ditches and 
their persisting in the present-day drainage system indicates consid- 
erable medieval investment in reclaiming this wetland. The two 
smaller enclosed compounds appear to centre around sandy areas 
within the larger enclosure (SI and S2 in Fig. 2). On the central 
platform within each compound, the K a measurements showed 
highly magnetic structures (Fig. 2b and d). At the easternmost site 

(51) , at least 18 features with a minimum diameter of 1 m, suggest 
the layout of a rectangular structure some 20 m by 30 m square 
(Fig. 2d). A central structure can be observed at the western site 

(52) (Fig. 2b). Here, the elevated K a of the ditches indicates their 
infilling with brick and other magnetic debris, whereas a trail of high 
K a values leading north suggests displacing rubble towards the sand 
ridge. 

To interpret and complement the geophysical data excavation 
trenches were laid out across the ditches and the central structure 
of each site (Supplementary Fig. 1-2). At both sites brick structures 
were found, as indicated by the K a measurements (Fig. 2b-e). These 
were identified as block foundations used to support a larger, most 
likely wooden, structure and to keep it away from the saturated soil. 
While a functional interpretation of S2 has been hampered by the 
poor conservation of the excavated building remains, the central 
structure at SI has been interpreted as a large monastic barn. The 
morphology of the sandy outcrop, combined with elevated ground 
water levels caused instability of the foundations leading to the col- 
lapse of the building. The excavated material culture and construc- 
tion materials indicate that both sites were occupied between the 13 th 
and the early 14 th centuries. This links their abandonment to a period 
of documented increased flooding in the region 21 . 

Based on the a a and excavation data, we modelled the relative 
medieval topography of the site 22 ' 23 . Although the heterogeneity of 
the ditch infillings, the peat layers and the targeted sandy layers 
limited the precision of this modelling procedure, an accurate model 
of the medieval surface was obtained. At 84 validation locations in 
the excavation trench at S2 a RMSEE of 0.37 cm and a Spearman 
correlation coefficient of 0.80 was obtained between the observed 
(Fig. 3a) and modelled (Fig. 3b) depths to the archaeological layers. 
Although the selection of the validation observations was not 
statistically random, the results confirmed the relative model corres- 
pondence to the medieval topography. A comparison of the mod- 
elled data to the observed surface along the validation transect 
(Fig. 3), shows the smoothing effect on the modelled data, but con- 
firms the accurate representation of the targeted archaeological 
layers. 

By integrating the K a data into the model, a comprehensive recon- 
struction of the reclaimed landscape was obtained (Fig. 4). Both SI 
and S2 are located on prominently higher and sandier locations, 
which are connected to the bordering sand ridge in the north. 
However, while the barn (at SI) was built on an existing sandy 
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Figure 2 | Electromagnetic induction data and excavated structures, (a), 1 m HCP a a data of the study area, plotted on an orthophotograph of the study 
area 33 , with indication of a man-made ditch (1), a sandy outcrop (2) and an area with accumulated peat and lacustrine deposits (3). The dashed lines 
indicate both moated sites SI and S2. (b), 2 m HCP K a data of the platform of S2, showing the central building and the magnetic disturbances caused by 
rubble in the ditches and leading towards the north (Supplementary Fig. 2). (c, e), Orthophotos of the brick foundations as revealed by the excavations at 
S2 (c) and SI (e). (d), 2 m HCP K a data showing traces of 18 brick block foundations on the platform in SI. Maps composed in ESRI ArcMAP 9.3.1. 



outcrop, excavations showed that S2 was artificially raised (Supple- 
mentary Figure 2d) by removing sand from the ditches and the adja- 
cent sand ridge. The design of both sites confirms a link to the nearby 
abbey buildings as their topographic position was chosen or altered 
in order to allow easy access to the sand ridge. The ditch enclosing the 
area was largest east of SI with a width of up to 10 m and a depth of 
1.5 m below the surface (Supplementary Figure Id). Here, a larger 
drainage capacity was needed, as the ditch was part of the complex- 
wide drainage system. This large ditch also expressed a separation 
between the inside and outside worlds; a symbolic reminder of the 
separation between religious and secular life 24 . The ditches circling S2 



and the remainder of SI rarely surpassed a depth of 1 m and had a 
width ranging from 2 m to 8 m (Supplementary Fig. Id and 2e). 
These wide ditches did not only supply additional drainage for the 
enclosed zones, but also formed a physical boundary that embodied 
aspects of identity and status 25 . Individual compounds, such as SI 
and S2, with an (artificially) raised platform and enclosed by ditches, 
can be defined as late-medieval moated sites 25 . In addition, the sur- 
face model shows dykes neighbouring the ditches that add to the 
embankment and visibility of both moated sites. Remnants of these 
earthworks were attested in the excavation trenches, either still in situ 
or thrown into the moat ditches (Supplementary Figure 2e). 
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Figure 3 | Medieval terrain model evaluation, (a), Geometrically correct visualisation of the western profile of the northern moat ditch in S2 
(Supplementary Fig. 2b). 84 calibration points, 0.2 m apart, along the targeted archaeological surface (orange line), were used to validate the surface 
model, (b), modelled medieval surface . Note the difference in y-scale between (a) and (b), which indicates the smaller standard deviation of the modelled 
depths (standard deviation = 0.10 m) compared to the observed surface levels (standard deviation = 0.43 m). 



Discussion 

Supported only by very limited invasive research, this reconstruction 
gives a comprehensive and unique insight into a designed medieval 
environment and shows, despite its transient nature, the amount of 
effort that was put into it. Our results also underline the discrepancies 
that can exist between historical information and uncovered archae- 
ological realities that can contribute to understanding the response of 
past societies to social and environmental changes. The 3-D mapping 
of multiple soil properties, combined with limited and directed 
invasive research, provides a broad foundation for further geoarch- 
aeological research. In the future, this methodology could be opti- 
mized by integrating additional physical soil properties, such as 
dielectric permittivity, with the a a and K a data. With the incorpora- 
tion of the presented approach into different studies and environ- 
ments, our knowledge of past human-landscape interactions can be 
significantly improved. 

Methods 

EMI survey. We used a motorized setup of a Dualem-21S EMI sensor 12 to survey the 
study area (Fig. 1). This EMI sensor simultaneously measures a a and K a with four 
coil configurations, i.e. with receiver coils at 1 m, 1.1 m, 2 and 2.1 m from the 
transmitting coil, in two orientations, i.e. horizontal coplanar (HCP) and 



perpendicular (PRP) orientation, resulting in 1 m HCP, 2 m HCP, 1.1m PRP and 
2.1 m PRP. The survey was conducted by driving across the study area along parallel 
lines, 0.75 m apart whereby measurements were taken every 0.25 m. During 
surveying, soil temperature was recorded at a depth of 0.3 m to correct the a a data for 
temperature variations 26 between survey days. The resulting measurements were drift 
corrected 12 . These data were then interpolated toaO.l m by 0.1 m grid using ordinary 
kriging 27 with the software Surfer v. 10 (Golden Software). All a a datasets were used 
but only the K a measurements from the HCP configurations were considered, as the 
K a data of the PRP configurations had a poor signal to noise ratio. 

Excavations. In August 201 1 and August 2012, two excavation trenches were laid out 
across the most apparent a a and K a anomalies (Supplementary Figures 1,2). All 
profiles and unearthed brick structures were recorded using computer vision-based 
3-D registration 28 . 

Depth modelling. The four a a measurements, representative for soil volumes down 
to 0.5 m, 1.0 m, 1.5 m and 3.2 m below the surface, were combined to model the 
depth to the soil layers below the more conductive plough layer and the peaty to clayey 
infilling of the ditches (Fig. 3). Based on the excavation data, the morphology of the 
targeted soil layers was considered representative for the medieval topography. We 
followed the procedure from Saey et al. 22 and De Smedt et al. 29 , based on the depth 
response functions of the Dualem-21S coil configurations as described by Wait 30 and 
McNeill 31 . Model calibration 23 was performed, using the observed depth to these 
layers obtained at 34 locations in the excavation trench at SI. 




Figure 4 | Model of the medieval landscape. Three-dimensional model of the depth to the medieval surface and schematic representation of the 
barn and central building remains at SI and S2, derived from the K a data (Fig. 2). The topographical model shows the sandy outcrop bearing structure SI 
( 1 ) and the artificially raised terrain forming the platform at S2 (2). At both sites, a pathway to the coversand ridge can be seen (3). In the north of the site 
the edge of the coversand ridge itself is visible (4). Moat embankments are present around both SI and S2 (5) and a large enclosure embankment 
(6) surrounds the entire area. Present-day parcel ditches separate the different agricultural fields in the area (7). 
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